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Abstract 

Experimental tests of Newton law put stringent constraints on potential deviations from standard 
theory with ranges from the millimeter to the size of planetary orbits. Windows however remain 
open for short range deviations, below the millimeter, as well as long range ones, of the order of 
qq ■ or larger than the size of the solar system. We discuss here the relation between long range tests 

m '. 

. of the Newton law and the anomaly recorded on Pioneer 10/11 probes. 
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I. INTRODUCTION 



The validity of Newton force law has been tested at distances ranging from the millimeter 

n q 

in laboratory experiments |1] to the size of planetary orbits. [2] But windows remain open for 
violations of the inverse square law at short ranges, below the millimeter, as well as long 
ones, of the order of or larger than the size of the solar system. This is also true for tests of 
general relativity which tightly constrain potential violations in the solar system [3| but let 
open space for deviations at small or large scales. 

The accuracy of short range tests has recently shown impressive progress, for gravity 
experiments pushed to smaller distances 
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ias recently shown impressive progress, for gravity 
> 4], 3 3 as wen as f° r Casimir force experiments. 3, 



Il2j In both cases, the agreement between theory and experiment is good. For 
Casimir experiments, it reaches an accuracy near the 1% level, after having accounted for 
the effects of imperfect reflection of the metallic mirrors used in the experiments. Q| 
The agreement can then be translated into constraints on potential violations of Newton 
force law at ranges from nanometer to millimeter, [to! llfij ] 

On the other hand, long range tests of the Newton law are performed by monitoring the 
motions of planets or probes in the solar system. The tests bearing on the third Kepler law 



or the precession of the perihelion of planets [l7( confirm the validity of general relativity. [181] 
The accuracy is especially good for ranges of the order of the Earth- Moon or Sun- 
Mars distances (see for example Refs. 20, 21, 22, 23). However, the Doppler data recorded 
on the Pioneer 10/11 probes show an anomaly when compared with calculations based on 
general relativity. |2J| The anomaly may be represented as an anomalous acceleration directed 
towards the Sun with a roughly constant amplitude. 25] It has not been explained to date 
though a number of mechanisms have been considered to this aim (see the discussions and 
references in Refs. 26, 27, 28, 29). 

At even larger scales, the rotation curves of galaxies show a conflict with general relativity 
as soon as the source of gravity is identified with the matter detected by electromagnetic 
means. Due to the excellent agreement of gravity tests with general relativity, this anomaly 
is commonly accounted for by introducing "dark matter" components designed to fit the ro- 
tation curves while keeping the gravity laws untouched. As long as dark matter components 
are not detected by other means, the anomaly can as well be ascribed to modifications of 
gravity at galactic scales. 

Q, 0, S H Q 

Similar statements apply for the observation 
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of an accelerated expansion through the relation between redshifts and luminosities for su- 
pernovae. This observation can be given alternative descriptions in terms either of "dark 
energy" or of modified gravity at cosmic scales. 35( Note that modifications of gravity are 



expected to be produced by vacuum induced effects |36j] or by effective gravity at low energy 
deduced from unification models. (3?], I3| An important requirement to be met by any such 
modification of gravity is that it is compatible with observations on galactic or cosmic scales 
while still matching the strict bounds set by gravity tests in the solar system. 

The Pioneer anomaly may be a central piece of information in this context by pointing 
at some anomalous behaviour of gravity at scales of the order of the size of the solar system. 
In the following, we focus our attention on the key question of compatibility of the Pioneer 
anomaly with other gravity tests. After having briefly recalled the observations, we consider 
the idea that the anomaly could be explained simply from a long-range deviation from the 
Newton potential, for example with a Yukawa form. We show that this explanation cannot 
be upheld against the data known for planetary tests of Newton law. More precisely, if the 
anomalous acceleration is ascribed to a Yukawa deviation from Newton law, the deviation 
is so large that it cannot remain unnoticed on the motions of outer planets, primarily Mars. 
This conclusion was already drawn [39] but it is written here under a form allowing us to put 
emphasis on the challenge raised by the incompatibility. The present paper also prepares 
discussions of a modification of Einstein theory of gravity to be presented elsewhere. |40| 



II. PIONEER ANOMALY 



The anomaly is recorded on radio tracking data from the Pioneer 10/11 probes during 
their travel to the outer parts of the solar system. At distances r from the Sun between 20 
and 70 astronomical units (AU), the Doppler data have shown a deviation from calculations 
based on general relativity. The anomaly is observed as a linear variation with time of the 
Doppler residuals, that is the differences of the observed Doppler velocity from the modelled 
Doppler velocity (see Fig. 8 of Ref. 25). It may be represented as an anomalous acceleration 
directed towards the Sun with a roughly constant amplitude on the range of distances over 
which it has been detected 

a P ~ 8 x 10" 10 ms- 2 , 20 AU < r < 70 AU. (1) 
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The anomaly can also be represented as a clock acceleration with the striking feature that 
its value ap/c ~ 3 x lCT 18 ss~ 2 is nearly equal to the Hubble frequency with a value H ~ 
80kms _1 Mpc _1 for the Hubble constant. 25 1 

Though a nntnber of —n. have been considered to this ai m ,0 H Q no 
satisfactory explanation of the anomalous signal has been found to date. Potential systematic 
effects do not seem to be able to reach the magnitude of the observed anomaly. Present 
knowledge of the outer part of the solar system does apparently preclude interpretations in 
terms of gravity or drag effects of ordinary matter. The inability of explaining the anomaly 
with conventional physics has given rise to a growing number of new theoretical explanations. 
It has also motivated proposals for new missions designed to study the anomaly and try to 
understand its origin (see the references in Refs. 26, 27, 28, 29). 

The importance of the Pioneer anomaly for fundamental physics and space navigation 
certainly justifies it to be submitted to further scrutiny. On the theoretical side, the incom- 
patibility of the Pioneer anomaly with other gravity tests appears to be a key question. We 
now discuss this question by considering the possibility that the anomaly could be explained 
from a long-range deviation from the Newton potential. To this aim, we use the common 
model of a Newton potential modified through the addition of a Yukawa perturbation. Q] 



III. YUKAWA MODIFICATION OF NEWTON LAW 

Throughout the paper, the potential energy V is written as the product of the mass 
energy of the probe by a dimensionless potential V (r) = mc 2 $ (r). The latter is the sum of 
the Newton potential and of a Yukawa correction 

•M = -%£(i + <w-i), (2) 

Goo is the effective Newton constant at large distances and M the mass of the Sun; A is the 
range of the Yukawa potential and its amplitude measured with respect to Gqo. The 
influence of the Yukawa perturbation disappears at the long distance limit r ^> A but it is 
significant otherwise. In particular, it gives rise to a potential linear in the distance r in the 
domain r -C A and, then, to a constant anomalous acceleration (see the next section). In 
the general case, the acceleration may be written as 



The anomalous acceleration contained in this expression could account for the Pioneer 
anomaly, if the range A is larger than 70 AU. But the value of the correction would thus 
be too large to remain unnoticed on planetary tests. In the following, we will put this 
incompatibility under a more precise form. 

Before going along this discussion, we rewrite the modified potential (J2J) in terms of a run- 
ning gravitational constant. To this aim, we write the Laplacian of its two components which 
obey respectively a Poisson equation (massless field propagation) and a Yukawa equation 
(massive field propagation) 

A (-1) = 4vr5(x) , (A - 1) (-Li) = ^(x). (4) 

The right hand sides represent point sources with 5(x) the Dirac distribution for 3- 
dimensional space position x. This equation leads to an expression of $ in Fourier space as 
a function of the spatial part k of the wavevector 

This expression may equivalently be written in terms of a running constant which replaces 
the Newton constant in the Poisson law 

G flcl M ~ f k 2 \ 

-k 2 $ [k] = 4vr^-L , G [k] = G 00 (l + a 00] -^ x J. (6) 

The running constant G [k] has been defined as the function of momentum and it goes from 
the value G m for small values of kA to the value (1 + a^) for large values of kA. Note 
that G(r), which can be defined from G [k] through an inverse Fourier transform, obeys 
A$ (r) = 4vr%^ (but not $ (r) = -^g^). 

IV. LONG RANGE LIMIT 

In the following, we will consider with a particular attention the long range limit, that is 
the domain r C A which could provide us with an explanation of the nearly constant value 
of the anomalous acceleration. To this purpose, we introduce an expansion of the modified 
potential (J2J) better adapted to this domain 

-r / \ GM GMa . _r S 

* (r) = --^ + MKr) , M>(r) = — j-fl-e *). (7) 



The gravitational constant and Yukawa amplitude have been redefined 

ot 

G = (1 + "oo) , Gat = GooOoo , a = t— — — • (8) 

1 + CKoo 

Note that G is the effective gravitational constant in the domain r C A. A power expansion 
of (JZJ) in terms of r/A leads to 

GMa 



<5$ (r) 



/ 2 

2A + °U 



(9) 



c 2 A 

The first term is constant and has no effect while the second one produces a constant 
anomalous acceleration. 

In general, the anomalous acceleration can be defined from (J3J) as 

a , \ GM . . . . . . . . GMa r _ t / r\i ... 

A(r) = — — + <L4(r) , (r) = — — |l - e a (^1 + -J j . (10) 
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It is then expanded in the domain r < A as 

. , , , GMa 

As already stated, the anomalous acceleration 5 A is essentially a constant in this domain. 
Should we identify it with the anomalous acceleration observed on the Pioneer probes, we 
would obtain an interpretation of the Pioneer anomaly. The sign convention is such that 
8 A < is needed to fit the observed anomaly. Note that the constant value of 8 A depends 
only on the combination a/ A 2 of the two parameters entering the Yukawa potential. Hence, 
there is a degeneracy in the choice of these two parameters which can be characterized by 
the relation 

^ = -^ , A = ^ 6300AU ~ 10 15 m. (12) 

A z A z V 2ap 

A is the distance from the Sun at which the Newton acceleration is twice the Pioneer 
anomalous acceleration. When plotted on a (A, a) diagram with log-log coordinates, the 
relation (JT2~j) defines a straight line with a slope 2. Note that only points with A > 70 AU fit 
a constant anomalous acceleration on the range of distances probed by Pioneer 10/11. Now 
the half-line thus defined ^70 AU < A, a = —j^J is excluded by existing planetary tests of 
Newton law. [2] The point is illustrated on Fig. 1 which shows recently updated constraints 
drawn from the analysis of the motions of planets and probes in the solar system: [18] the 
values (A, a) fitting the anomaly are clearly in the forbidden domain. This conflict will be 
written directly in terms of the anomalous acceleration 8 A in the next sections. 
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FIG. 1: Constraints on a Yukawa modification of Newton law. The colored domains represent, on 
a log-log plot, the excluded values for amplitude a as a function of the range A (measured in m). 
Courtesy of the authors of Ref. 18. 

The discussion can also be translated into terms bearing on the running constant. To this 
aim, we rewrite © under a form better adapted to the domain kA ^> 1 (which corresponds 
to r < A) 

G[k] = G + 5G[k] , ^M = -k^T ( 13 ) 
We then expand the anomalous part in this domain 



SG [kl 



Ga 

~k2p 



1 + 



k 2 A 2 



(14) 



As the constant anomalous acceleration S A in (fTTj) with which it is directly associated, the 
anomalous term 5G [k] is proportional to the combination a/A 2 . Note that Eqs. Q, (jlljl 
and (JHJ), valid in the domain r <C A (or equivalently kA 3> 1), are sufficient for the purpose 
of our discussions. At the same time, they have been derived from Eqs. (J2J), (jUJ) and (JHJ) 
which are better behaved at the limit of large distances r ^> A (or equivalently kA <C 1). 



V. PLANETARY TESTS 



Tests of the Newton law bearing either on the third Kepler law or on the precession of 
the perihelion of planets are known to confirm the validity of general relativity with a good 
accuracy. We now make the significance of this statement more explicit by writing directly 
in terms of the anomalous acceleration 5 A the constraints drawn from these tests. Since we 
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study gravity in the outer solar system, we use the Newton theory with the Sun described as 
a motionless point source. The motion of the probe mass takes place in the central potential 
(J2J) and we can write the conservation of energy E and angular momentum J 

E V r + V l , * ( \ 3 
2 = o 2 + $ r , — = ™* 

' u * = r dF' (15) 

r is the distance from the Sun, i the time coordinate and (p the azimutal angle; the trajectory 
has been assumed to stay in the plane 9 — | . When eliminating time, we obtain the equation 
of motion as 

d 2 u m 2 2 . , 1 

— — + u = — — rA(r) , u=-. (16) 
etyr J z r 

Note that r 2 A (r) can be replaced by c 2 |^ when $ is the potential expressed in terms of the 
variable u. 

Using the standard Newton law for the potential (a = 0), the right hand side in (jlfi)) is 
a constant and the Kepler ellipse is recovered 

l + ecosy? ( 2X J 2 GMm 

M= ^^ ' P = a (1 — e ) — -QpJ^2 • (17) 

The orbital frequency cu = with T the orbital period, is given by the third Kepler law 
LJ 2 a 3 = GM or, equivalently, by an expression directly drawn from (jlfij) 

m2 r 9 .i GMm 2 , , 

M- = -ja t^L = -72-- ( 18 ) 

Here [u] st and [r 2 ^4] st are the values obtained for u and r 2 A in standard theory. 

With the modified Newton law (a ^ 0), there is a correction to the third Kepler law (|18jl. 
If we evaluate it on circular orbits (e = 0), we obtain in a linear approximation with respect 
to the small perturbation 

u r 2 A r 2 5A(r) 
Consider now a planet, say Mars, for which the elements a or p describing the orbit have 



been measurec 
orbital period. 



through optical means, that is independently of the measurement of the 
171 ] The result obtained for Mars can then be compared with the reference 
value fixed by the orbit of Earth by forming the ratio WMars/wEarth whose difference to unity 



now depends on the variation of ()19|) from Earth to Mars. Defining the relative accuracy 

8 



FIG. 2: Sensitivity of a Kepler test (left hand part) and precession test (right hand part). The full 
lines delineate the exclusion domain (above the curve) for the value of the reduced amplitude a/e 
(definitions of e given in the text) as a function of A (measured in AU). The logarithms go from -2 
to +6 on the horizontal scale and from -1 to +9 on the vertical scale with the zeros indicated by 
thin lines. 

e = \Sa/a\ on the element a of Mars as it is measured in AU, we deduce that an anomalous 
acceleration could be noticed under the condition 

(20) 





~a 2 5A (a)' 




~a 2 5A (a)' 
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Using (JHJ), this condition is read as a minimum value for the amplitude 

1 



a 

— > 

e 



(21) 



[(i + ^^Mars-^ + ^^Earthl' 

The curve showing the frontier of the domain is plotted on Fig. 2 with log 10 — log 10 scales. 
The range A as well as the radii are measured in AU (a = 1 for Earth and 1.5 for Mars). 
The long range asymptote on the curve corresponds to a fixed value for a/ A 2 and, therefore, 
for the anomalous acceleration S A. This point will be discussed in more detail below. 

Precise planetary tests are also performed by following the precession of the perihelion 
of planets. |2j, ll2| For simplicity, we evaluate it for an orbit with a low eccentricity e < 1 
in a linear approximation with respect to the small perturbation. In this case, the variable 
u undergoes a small sinusoidal variation around a constant 1/p ~ 1/a (see Eq. ITTj) . The 
equation of motion (jlti)) can thus be replaced by an approximation corresponding to the 
linearization of the associated variation of r 2 A (r) 



The parameter 1/p has been replaced by a modified value 1/P but this does not matter 
for our present purpose. What is important for the evaluation of the precession is the 
modification of the coefficient in front of u. It implies that the perihelion (maximum value of 
u) is recovered when tp has ran over 2ir (1 + 8k) rather than the standard value 2tt. Defining 
the relative accuracy e = \8ip /27r\ for a test of the precession, we obtain the following relation 
for the anomaly to be detectable 

(23) 

Mars 

The frontier of the exclusion domain is drawn as the right hand plot on Fig. 2. It has roughly 
the same shape as the left hand plot and, in particular, the long range asymptote on the 
curve corresponds again to a fixed value for a/X 2 . 



e < a 



2A 2 



VI. DISCUSSION 



As already stated, the Yukawa correction of the Newton law would produce a constant 
anomalous acceleration over the range of distances probed by Pioneer 10/11 provided the 
Yukawa range is large enough A > 70 AU. This means that the tests performed on planets 
or probes in the solar system correspond to the limit of long ranges r C A. As shown in 
the preceding section, this also means that they test a single combination a/X 2 of the two 
parameters entering the expression of the Yukawa correction. This entails that it is possible 
to write the constraints directly in terms of this anomalous acceleration 8 A. 

This fact is especially clear for the Kepler test since the condition (|2T)|) can be rewritten 

8 A ((ZMars) 8 A (a E arth) 



£ < 



(24) 



A (ctMars) A (dEarth) 

Now 8 A (a) has essentially the same value at a = a Ma _ rs and a = aEarth as well as at any 
distance smaller than the Yukawa range. If we simply denote this constant anomalous 
acceleration as 8 A, we deduce that a test with a relative accuracy e is immediately translated 
into a bound on 8 A 



\8A\ < e 



^Earth^Mars 



(25) 



^Earth — ^Mars 

Considering 39] that the distance to Mars has been tested with an accuracy of the order of 

10 -10 . Inserting the values of ^Mars ~ 2.6 x 10~ 3 ms -2 and ^Earth ~ 



12m, we obtain e 
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5. 9 x 10~ 3 ms -2 , we deduce that 5 A should remain smaller than 5 x 10 -13 ms -2 . This is 
certainly much smaller than the value needed to explain the Pioneer anomaly. 

For the perihelion test, the accuracy (|23|) can be rewritten in the long range limit r <C A 

as 

Am A A i n. t I 

(26) 





5ip 




5 A (a M ars) 






2tt 
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v OMars) 




It follows that a test with an accuracy \5ip\ 


is again translated into a bound on 5 A 




5A 


< 


2tt 


| ^Mars 





(27) 

Both results (J25)) and (|27|) constrain the value of \SA\ that is also the combination a/A 2 of 
the two Yukawa parameters. One can therefore extract the best bound on \5A\ from the 
long range asymptote on the diagram (A, a) of Fig. 1. This diagram collects the most recent 
planetary data 



181 ] and it leads to the following bound 
<L4 < 3 x 1(T 13 ' 



ms 



(28) 



We can now sum up the present paper as discarding the possibility that the Pioneer 
anomaly could be explained from a long-range Yukawa deviation from the Newton potential. 
This conclusion 39] has been written directly in terms of the anomalous acceleration 5 A 
which appears in such a modification of Newton law. The value needed to fit the Pioneer 
anomaly is in fact more than 1000 times too large to remain unnoticed on tests of the 
Kepler law or of the precession of perihelions. The discrepancy illustrates the challenge to 
be met when trying to analyze the Pioneer anomaly in the same framework as other tests 
of gravity. j3] 
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